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Abstract. Collision processes between two gold clusters are investigated using classical molecular dynam- 
ics in combination with embedded atom (EA) potentials, after checking the reliability of EA results by 
contrasting them with first principles calculations. The Au projectiles considered are both single atoms 
(N=l) and clusters of N=2, 12, 13 and 14 atoms. The targets contain N= 12, 13 and 14 gold atoms. The 
initial projectile energy E is in the range < E < 1.5 eV/atom. The results of the collision processes are 
described and analyzed in detail. 

PACS. 36.40.Qv Stability and fragmentation of clusters - 
36. 40. Mr Spectroscopy and geometrical structure of clusters - 

61.46.+W Nanoscale materials: clusters, nanoparticles, nanotubes, and nanocrystals - 

73.22.-f Electronic structure of nanoscale materials: clusters, nanoparticles, nanotubes, and nanocrystals - 

82.30.Nr Association, addition, insertion, cluster formation 

1 Introduction has also come to the forefront of technology On the 

theoretical side ab initio procedures are now capable of 

The study of nanostructures has recently attracted wide- providing incisive insights into the properties of these sys- 

spread interest among theoretical and experimental physi- terns. In addition, novel and sophisticated nanostructure 

cists and chemists, and because of its many applications fabrication, manipulation and measurement techniques have 

o j ~Z ~, 7 ~7~ hi- i tv • given impetus to experiment, and reliability to a large 
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stimulated activity in the nanostructure field since they proven to be quite reliable pC[| . Classical molecular dy- 

can be considered as building blocks of novel nanostruc- namics (MD) is a valid option in this case, since it allows 

tured materials and devices ^,0,EolO]. to considerably reduce the computation time (relative to 

In particular, metallic clusters provide an interesting MD ab lmUo calculations), and/or to increase substan- 

subject of study for at least two reasons: i) clusters con- tiall y the number of particles that can be handled. Ob- 

stitute intermediate systems between isolated atoms and viousl y there is a P rice to P a ^ as examples of this cost 

molecules, on the one extreme, and bulk solids on the we mention that the method is at best semi-classical, that 

other {i.e. they constitute genuine mesoscopic systems); the detailed electronic structure is ignored, and that phe- 

and ii) often they exhibit an interesting phenomenology of nomcnological potentials, adjusted to fit bulk properties, 

their own. Gold clusters have received widespread atten- are used. In spite of these shortcomings the MD technique 

tion during the last two decades, both experimentally §, can reliably be used to compute the properties dominated 

and theoretically 0,[||l|,|l|,|l6|,|l7[|l|,|l|,|o[ b y th e ionic contribution, which is the case for the phc- 



2l] p2||23| , |24| ^5[| nomena treated in this paper. This is especially true when 

MD is checked, in some physical limit, against first prin- 



ciple calculations. An alternative, but equally convenient 
procedure, is to close in on a solution, for example for 
geometrical optimization via MD, to be followed by first 
principle calculations. The synergy between ab initio and 
MD thus allows to significantly reduce the resources that 
are required and to expand the set of problems amenable 
to treatment. 



In principle one does expect ab initio procedures to be 
the definitive tool to handle this type of systems; how- 
ever, it is not always feasible (or at least practical) to 
implement ab initio calculations. In fact, since we are in- 
terested in structures with a fairly large number of atoms, 
or arranged in several different interacting nanostructures, 
ab initio computations first become very time consum- 
ing and, in the end, impractical. Much the same happens 
when attempting to obtain a detailed description of the More than a decade ag0 the bombardment with gold 
long time evolution of small systems, or to describe their clusters f metallic surfaces was investigated experimen- 
properties when subject to a large variety of external con- tally and theoretically @||,||. Moreover, the dra- 
ditions. matic energy accommodation that occurs in cluster-cluster 

Here we intend to develop an adequate description of collisions, which is crucial to understand the growth mech- 
colliding gold clusters, a phenomenon which falls into one anism during the early stages of particle formation, was 
or more of the categories described in the preceding para- investigated around the same time by Blastein et al. [ j30| . 
graph. In this paper we employ mainly classical molecu- We focus our interest on the dynamics of gold cluster colli- 
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of small Au clusters, and the dynamics of atom-cluster 
and cluster-cluster collisions, for different values of the 
impact parameter b and as a function of center of mass 
energy E. 

This paper is organized as follows: after this Introduc- 
tion we describe, in Sec. |^, both the ab initio and EA meth- 
ods employed in the computations. In Sec. ^ we provide 
results, for a large variety of cases, of the implementation 
of the codes, for the cluster structures, their symmetries 
and the dynamics of the collision process. Finally, in Sec. [| 
we draw conclusions and close the paper. 

2 Simulation method 
2.1 Ab-initio method 

In order to asses the quality of the semi-empirical em- 
bedded atom (EA) procedure used in the context of our 
classical molecular dynamics simulations, we start con- 
trasting EA against first principles calculations. Thus, ab 
initio geometrical optimization was carried out within the 
Car-Parinello approach [pl[ , in the framework of the den- 
sity functional theory, using gradient corrections in the 
PBE implementation p^ |. Gradient corrected functionals 
have been adopted in recent theoretical studies of geo- 
metrical optimization of metallic clusters, mainly because 
they are more accurate than the local density functional, 
even though there still is some controversy in the literature 
on the ground state geometry of small Au clusters |l5|,j33"|, 
. The calculations were performed only at the /"-point 
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tentials [ ^5| . The wavefunctions were expanded in plane 
waves, with an energy cutoff of 60 Ry. We have explic- 
itly checked that, with this energy cutoff, the structural 
properties of our system are well converged. The box used 
in the calculations was always at least three times larger 
than the cluster diameter. The results obtained, as well 
as their EA counterparts, are given in Fig. |^ and in Ta- 
ble ^. The ab initio optimized geometries are illustrated in 
Fig. |]. The agreement between ab initio, EA and Wilson 
and Johnston |0J (WJ) procedures is quite satisfactory 
for large clusters (N<10) and provides a reasonable basis 
to trust the EA calculations that constitute the core of 
the present paper. 
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Fig. 1. Ab initio optimized geometries for 3, 4, 6, 8, 10 and 12 atom gold clusters. 
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2.2 Embedded atom method 
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The interatomic interaction between gold atoms is mod- The first issue we address is to check the reliability of 



eled using a semi-empirical (EA) potentials |36j|37|]. On our MD procedure as compared with alternative meth- 

the basis of these EA potentials we obtain the average ds. In Table § we compare the average binding energy 

binding energy per atom, E b , which is later minimized to E b of the lowest energy configurations, and the average 

yield the optimal cluster geometry. The latter is achieved nearest neighbor distance R, of clusters built with dif- 

using a Monte Carlo procedure, for which we adopted as ferent numbers of Au atoms. On the one hand we have 

starting configurations, for the different cluster sizes, the ab initio results obtained within the framework of density 

geometries found by Wilson and Johnston ||l| (WJ). Once functional theory, and on the other the empirical potential 

the optimal geometry is established several static proper- results of WJ Q (who used the Murrell-Mottram 

ties, like nearest neighbor distances and angles, and the potential) as well as the EA values that we obtained, 

average coordination number, are readily evaluated. R ig quite apparent that the EA estimates for Eb diffcr 

After the various different clusters are properly charac- considerably from the ones found by WJ, and also with 

terized we use classical molecular dynamics (MD) to sim- the more trustworthy ab initio results. However for small 

ulate the cluster-cluster scattering process. Many body clusters the geometrical parameters are not so satisfac- 

EA semi-empirical potentials are used throughout. To in- tory. This small N (where N is the number of atoms in the 

tegrate the equations of motion we implement the Verlet cluster) error margin is not unexpected, since the EA po- 

velocity algorithm, with a 1 femtosecond time step. Since tential has been adjusted to fit bulk properties and cannot 

the collision fragments heat up as a consequence of the be expected to fully succeed in systems where N is tiny, 

scattering process it is necessary to cool them down; to However, and also as expected, the situation improves as 

do so we simply rescale the temperature by 1 % every N increases, which is precisely what is borne out by Ta- 

1000 steps, during a total of 100 000 time steps, to reach ble g, where we display results for the range of clusters 

a final energy of 63 % of the initial one. Finally, the col- s i zes from N=3 to 12. The same trend is observed in Fig. §, 

lision fragments are carefully scrutinized to extract the where we plot E b versus N. We notice that ab initio cal- 

physical information we arc looking for. culations were performed for N = 4, 5 and 6 planar and 

three dimensional structures. It is apparent that as the 
cluster size increases, specially when the number of atoms 

3 Results and Discussion is l ar g er than 10, that both the ab initio and EA values 
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Fig. 2. Binding energies per atom Eb, obtained by WJ using 
the Murrell-Mottram potential, and by us on the basis of the 
EA potential and ab initio, as a function of the number of 
atoms in the cluster N. The ab initio results, for which we have 
considered several possible geometries, are detailed in Table ^. 

The average nearest neighbor distances R that we com- 
pute also exhibit larger errors than those of WJ. In spite 
of the fact that we employ the WJ cluster configurations 
as the starting point for our calculations, but in which we 
use a binding energy obtained from a different potential, 
we derive geometrical structures which differ from those of 
WJ. However, once again, increasing N yields compatible 
results. For example, for a 3-atom cluster the difference in 
nearest neighbor distances amounts to 15%, but already 
for a 6-atom cluster it reduces to only 8%. These aver- 
age distances are illustrated in Fig. pi where we plot R as 
a function of the number of atoms in the cluster, N. In 
contrast with Fig. ^ the plot of Fig. || is not smooth, but 
shows abrupt variations between two successive values of 
N. Despite this roughness the tendency of the WJ and our 
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Fig. 3. Average nearest neighbor distances R calculated using 
EA (dashed line) , and those obtained by WJ using the Murrel- 
Mottram potential (full line). 

our results are in good agreement with the ab initio ones 
obtained by Wang [^6j . 

3.2 Symmetries 

In addition to the binding energies and interatomic dis- 
tances the cluster symmetry is a relevant characteristic 
and, in the context of gold cluster topologies, the Jahn- 
Teller effect is also an important element. Ab initio calcu- 
lations predict a C2„ symmetry for AU3 and Au4, while EA 
yields C2 V and T>2d, and the Murrell-Mottram potential 
used by WJ yields D^h and T^ symmetries, respectively. 
These differences are quite apparent in the pair correlation 
function g(r) plotted in Fig. || . 

The second difference in the binding energy is defined 

by 
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Fig. 4. Pair correlation function g(r) for 20, 30 and 40 Au 
atom clusters as calculated by us, using EA potentials, and by 
WJ. 

and gives an indication of the stability of a cluster with 
respect to disproportionation Q, as well as its ionic hard- 
ness. A plot of A2Ef J (N) as a function of N is given in 
Fig. ||, where we observe a good agreement of our EA val- 



ues with those reported by Wilson et al. 21 , both for 
the position and magnitude of the hardness peaks. The 
maxima (minima) of A%E})(N) imply that there are val- 
ues for which it is more difficult (easier) to add an atom 
to the cluster. Moreover, as N ^> 1 the plot becomes quite 
smooth, which constitutes an indication that the cluster 
can incorporate an additional atom without major hin- 
drance. 



0.4 



0.2 



> 



-0.2 




Fig. 5. Second difference in the binding energy Zi2{£*(N)} as 
a function of cluster size N. 

3.3 Collisions between a single Au atom and a Au 
cluster 

Next we report the results of our simulations of the col- 
lisions between a single Au atom and a variable size Au 
cluster, for several values both of the initial per atom en- 
ergy E and of the impact parameter b. The precise details 
that describe the collision process are as follows: at time 
t = the centers of mass of the atom and the cluster 
are placed along the x-axis, ±10 A away from the origin, 
respectively. The atom is located a distance b away from 
the x-axis, along the y-axis on the xy plane. The principal 
symmetry axes of the various clusters are aligned perpen- 
dicular to the direction of motion, that is parallel to y. The 
impact parameter b is varied between (head-on collision) 
and 7 A; the latter corresponds to the distance where the 
interaction potential effectively vanishes, since the average 
radius of a cluster with 12 < N < 14 is less than 3 A and 
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varied between 0.1 and 1.5 eV per atom, in steps of 0.2 eV on the basis of the following symbols: single atoms, dimers 

(8 values in all). The maximum energy E = 1.5 eV per and trimers are represented by dots, two dots joined by 

atom corresponds to approximately one half of the cluster a line, and three dots that form a triangle, respectively. 

binding energy. After the collision takes place the result- When several of these collision fragments are generated 

ing fragments are stabilized, by gradual cooling through we denote their number by a factor in front of the corre- 

the rescaling of the internal velocities. Finally, we analyze sponding symbol. If the fragment contains four or more 

the data characterizing the collision fragments for several atoms the symbolic representation is a circle. Finally, if 

special cases. the number of atoms in the target is not altered after the 

,„ . . u , r • i collision, we stress this fact representing it by a square 

We consider two categories: i) the scattering of a single ° J n 

ii , • , , , tvt 10 10 jia A with the original number of target atoms in its interior, as 

gold atom against clusters with JN= 12, 13 and 14 atoms, & b ' 

which is dealt with here; and, ii) the scattering of a vari- l° n S as N> 4. 

able size projectile (N= 12, 13 y 14) on a variable size In Fig. || it is noticed that for low energies 

target with a similar number of Au atoms, which is pre- and small impact parameter (0 < b < 3 A), projectile 



scntcd in 3.4. Throughout we use the concepts of low and and target fuse into a 13 atom cluster. When the energy 
high energies, and small and large impact parameters. Low is increased to 0.9 < E < 1.5 eV, keeping the impact 
energies are defined to be in the range 0.1< E <0.7 eV parameter fixed, coalescence is observed for a few cases, 
and large within 0.7< E <1.5 eV. Similarly, small impact while more often 11 or 9 atom clusters, and one or two 
parameters cover the range < b < 3 A, and large is de- dimers, respectively, are generated. For larger impact pa- 
fined as 3 < b < 7 A. The upper bounds on E and 6, as rameters, 4 < b < 7, fusion is present in the low energy 
mentioned above, are related to the binding energy of the region, but for a few cases target and projectile remain 
cluster and the distance at which the projectile does not unaltered. Finally, for large b and E values we observe 
interact with the target, respectively. a few cases of coalescence, some 11 atom clusters plus a 

^ , ion f j.u • i r u dimcr and many instances (denoted by squares) in which 

One atom on 12 Because ol the rich variety of results, J v J H ' 

, . r , ,i • i , j. i ,i i projectile and target size do not change, 

and to facilitate their understanding by the reader, we F J ° ° 

have chosen to illustrate them by means of figures. In par- A particularly interesting scattering process, illustrated 

ticular Fig. ^ describes the one gold atom collision with a in Figs, [t], occurs for an energy of 4 eV and impact param- 

12 atom cluster. In this and ensuing figures the columns etcr 1.8 < b < 2.2 A. When the projectile approaches the 

correspond to different values of the energy per atom (in target (Fig. |?]a) it attracts the nearest (lowest in Fig. f?]b) 

eV), while the rows correspond to several different impact atom, but without removing it from the cluster. This gen- 
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Fig. 6. One Au on 12 atom cluster collision. The numbers on 
the left column denote the impact parameter b, measured in 
A, and the top row the average energy per atom E in eV. 

amplitude vibrations in the "lowest" atom, as the projec- 
tile leaves the scene (Fig. 0c). As a consequence of these 
large amplitude vibrations this particular atom does over- 
come the energy barrier and ends up at the center of the 
cluster (Fig. ^c). We have estimated a lower bound for 
this energy barrier of 0.03 cV. Quite remarkably this lower 
symmetry cluster has an energy slightly smaller than the 
fully symmetric 12 atom cluster we accepted above as the 
stable configuration. 



o 



(a) 




(b) 



O 



Fig. 7. Collision between a single Au atom and a 12 atom 
cluster. The single atom approaches the cluster from the left 
(a), and exchanges energy (b). Finally, in (c) one of the cluster 
atoms is displaced towards the target center, while the projec- 



In fact, the original 12(1 h) cluster has a binding energy 
Ei, ~ —3.03 eV, while the energy of the less symmetric 
configuration 12(C5„) equals Ef, = —3.09 eV, which is 
2% lower. Thus, we are faced with the question of why 
this lower energy configuration was not obtained in the 
minimization process we reported above. The explanation 
of this apparent contradiction is related to the fact that 
the energy barrier the atom has to overcome, to shift to 
the center of the cluster, is rather large and cannot be 
achieved in a minimization process that starts with an 
icosahedron and allows only small displacements from the 
original equilibrium positions. It is worth mentioning that 
this asymmetric structure was also asigned minimum en- 
ergy by Wilson and Johnston |2l). However, they argue 
that the bond compression of the 12(C5„) structure gen- 
erates a repulsion strong enough to destabilize it in favor 
of the icosahedron. 

One atom on 13 When a Au atom collides with the 
"magic number" N=13 cluster the results are not much 
different from the previous 1 on 12 case. However, there 
is a larger number of coalescence cases, many instances of 
dimcr formation, and also instances are observed where 
there is no change in the number of atoms of projectile 
and target. For example, when E = 0.9 eV, there is fusion 
for b < 2, dimer plus a 12 atom cluster generation for 
b — 3, fusion for b — 4, no change b — 5, and fusion for 
b = 6 and b = 7 A. 

One atom on 14 In this case the results fall into only 
two categories: either target and projectile fuse or they 



10 Jose Rogan et al.: Rearrangement 

For b > 5 A there is coalescence only for very small en- 
ergies E < 0.3, and no change in the number of atoms in 
projectile and target for larger energies. This small 6 large 
E behavior can be understood as follows: for low impact 
parameter the collision gives rise to violent cluster vibra- 
tions and deformations, which precludes the trapping of 
the projectile. Instead, larger b values induce less dras- 
tic cluster deformations and sufficient attraction to fuse 
projectile and target. 

One atom on 12: cluster rotation In all preceding 
cases the cluster symmetry axis was taken to be parallel 
to the y axis, that is, perpendicular to the initial projectile 
velocity. Now we align the cluster symmetry axis paral- 
lel to the y axis. A qualitative change is observed: single 
atoms as a result of the collision. In fact, for low ener- 
gies and small b there is fusion. However, for E > 0.7 eV 
several alternatives are observed: either coalescence, or 8, 
9 and 11 atom cluster formation accompanied by the cre- 
ation of dimers and single atoms. When b > 3 A we obtain 

11 atom clusters plus a dimer or a pair of isolated atoms. 
Finally, for the largest b — 7 A value no atomic reordering 
is observed. However, we notice that in general the overall 
structure of the results is equivalent to the perpendicularly 
oriented cluster impact discussed above, which allows us 
to concentrate on cluster collisions without paying much 
attention to their relative spatial orientation. 

3.4 Collisions between two Au clusters 

Now we turn to the problem of the collision of two clusters 
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One dimer on a 12 atom cluster 

Fusion is obtained for low energies (E < 0.7 eV) and 
practically all values of the impact parameter b. In many 
instances the end result is a dimer and a 12 atom cluster, 
which are the outcome after a complex dynamic inter- 
action that, finally, yields a reconstruction into the two 
original clusters. For E > 0.7 eV and small b we notice a 
diversity of results: single atoms plus clusters of 5, 6, 7, 8, 
10 and 11 atoms are obtained, which reflects the fact that 
the more complex the projectile the richer the variety of 
collision fragments. 

One 12 atom cluster on another 12 atom cluster 

Fusion is observed only for low energies (E = 0.1 eV 
per atom) over the whole range of impact parameters 
< b < 7 A. For large energies and small impact pa- 
rameter collisions, either large fragments plus a couple of 
dimers or trimers, or dimers and single atoms are gener- 
ated. For large E and large b, collisions without cluster 
size rearrangement are predominant. 
One 13 atom cluster on a 12 atom cluster Again 
there is coalescence for E ~ 0.1 eV and < b < 7 A, while 
for large E and large b collisions without rearrangement 
predominate. In the ~ b ~ 5 A and E > 0.5 eV region 
large and medium size fragments plus trimers, dimers and 
single atoms are produced. Total break up of the cluster is 
seen almost exclusively for small b and large E collisions. 
One 14 atom cluster on a 12 atom cluster Again 
here the results are quite similar to the previous ones. 
For small b there is fusion, while for large b and large 
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clusters. Moreover, for a large region of parameter space 
(0 < b ~ 5 A and 0.5 ~ E ~ 1.5 eV) a whole vari- 
ety of fragments does result: large fragments plus dimers, 
medium size fragments and finally, in for the largest ener- 
gies, total cluster breakup into small pieces. 
One 13 atom cluster on a 13 atom cluster 

The collision of two "magic number" clusters yields 
the rich variety of results illustrated in Fig ||. It is read- 
ily noticed that there are many notable exceptions to the 
general trends observed in the preceding cases, and which 
are only present for this particular case. This is specially 
noticeable in the upper right of the figure, which illus- 
trates the parameter values for which dimers, trimers and 
clusters of 6, 7, 8, 9, 11, 12, 14, 15, 16, 17 and 24 atoms 
are generated. However, again fusion is observed for low 
energy collisions (E ~ 0.1 eV) for all b values, as well as 
no rearrangement collisions for large impact parameters, 
which are the final outcome of a complex dynamic inter- 
action that in the end yields a reconstruction into two 
clusters equal to the original ones. 
One 14 atom cluster on a 13 atom cluster 

A low energy coalescence region, as well as a large im- 
pact parameter zone where after the collision projectile 
and target rearrange into their original structures, is again 
obtained. Also, and just as in previous cases, for large E 
and small b a diversity of fragments (dimers, trimers, and 
5, 6, 7, 8 and 9 atom clusters) do result. In addition there 
are several cases with a single Au atom exchange between 
projectile and cluster, such that the original (14, 13) pair is 
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Fig. 8. Collision between two 13 atom clusters. The numbers 
on the left column denote the impact parameter 6, measured 
in A, and the top row the average energy per atom E in eV. 

E region when three, instead of two, collision fragments 
are generated, always a dimer (and not a single atom) is 
created. This dimer originates in the 14 atom projectile, 
yielding in the end a 12 atom and a 13 atom cluster plus 
a dimer, as collision fragments. 

4 Conclusions 

The dynamics of the collision process of gold clusters has 
been investigated by means of classical molecular dynam- 
ics in combination with embedded atom (EA) potentials. 
First, the reliability of the EA potentials was confirmed 
by comparison with ab initio values, finding that EA is 
in good agreement for the cluster sizes we considered, an 
agreement which improves as N increases. Next, struc- 
tural characteristics and the symmetry of the various Au 
clusters were obtained and contrasted with published re- 
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Several type collisions were investigated, finding re- 
gions of coalescence, fragmentation and scattering. Which 
of these outcomes actually occurs depends mainly on the 
values of the projectile energy E and the impact param- 
eter b. Coalescence is dominant for low energies (E < 
0.7 eV) and small impact parameters (b < 3 A). Sim- 
ple scattering, with no change in the size and structure of 
the colliding clusters, prevails for large E and b values. For 
large energies and small impact parameters fragmentation 
and scattering are generally the case. For large E and large 
b scattering is the most probable outcome. When the clus- 
ter does breaks up the main collision products, apart from 
large fragments, are dimers. On the other hand cluster co- 
alescence provides a viable mechanism to generate larger 
cluster sizes. It is also of interest that the collisions them- 
selves turn out to be rather insensitive to the relative ori- 
entation of the projectile and target main symmetry axes, 
and that cluster collisions can generate metastable struc- 
tures, which are usually not accesible due to the existence 
of a potential barrier. 



collisions between gold clusters 

References 

1. W.A. de Heer, Rev. Mod. Phys. 65, 611 (1993). 

2. C.L. Cleveland, U. Landman, T.G. Schaaff, M.N. 
Shafigullin, P.W. Stephens and R.L. Whetten, Phys. Rev. 
Lett. 79, 1873 (1997). 

3. T.G. Schaaff, W.G. Cullen, P.N. First, I. Vezmar, R.L. 
Whetten, C. Gutierrez-Wing, J. Ascencio and M.J. Jose- 
Yacaman, J. Phys. Chem. 101, 7885 (1997). 

4. K.J. Taylor, C.L. Pettiette-Hall, O. Cheshnovsky and R.E. 
Smalley, J. Chem. Phys. 96, 3319 (1992). 

5. K. Koga, H. Takeo, T. Ikeda and K.I. Ohshima, Phys. Rev. 
B57, 4053 (1998). 

6. V.A. Spasov, Y. Shi, K.M. Ervin, Chem. Phys. 262, 75 
(2000). 

7. B. Palpant, B. Prevel, J. Lerme, E. Cottancin, M. Pellarin, 
M. Treilleux, A. Perez, J.L. Vialle and M Broyer, Phys. 
Rev. B57, 1963 (1998). 

8. R.L. Whetten, J.T. Khoury, M.M. Alvarez, S. Murthy, I. 
Vez-mar, Z.L. Wang, P.W. Stephens, C.L. Cleveland, W. 
D. Luedtke and U. Landman, Adv. Mater. 5, 8 (1996). 

9. R.P. Andres, T. Bein, M. Dorogi, S. Feng, J.I. Hender- 
son, CP. Kubiank, W. Mahoney, R. G. Osifchin and R. 
Reifenberger, Science 272, 1323 (1996). 

10. C.A. Mirkin, R.L. Letsinger, R.C. Mucic and J. Storhoff, 
Nature 382, 607 (1996). 

11. A.P Alivisatos, K.P. Johnsson, X. Peng, T.E. Wilson, C.J. 
Loweth, M.P. Bruchez and PG. Schultz, Nature 382, 609 
(1996). 

12. J.J. Zhao, X.S. Chen, G.H. Wang, Phys. Lett. A189, 223 
(1994). 

13. H. Handschuh, G. Gantefor, P. S. Becthold and W. Eber- 



Jose Rogan et al.: Rearrangement 

14. I.L. Garzon and A. Posada- Amarilla, Phys. Rev. B54, 
11796 (1996). 

15. O.D. Haberlen, S.C. Chung, M. Stener and N. R6sch, J. 
Chem. Phys. 106, 5189 (1997). 

16. J. Wang, G. Wang and J. Zhao, Phys. Rev. B66, 
manuscript #BMA8242, to be published Aug. 15, 2002. 

17. I.L. Garzon, K. Michaelian, M.R. Beltran, A. Posada- 
Amarilla, P. Ordejon, E. Artacho, D. Sanchez-Portal and 
J.M. Soler, Phys Rev Lett. 81, 1600 (1998). 

18. J.L. BelBruno, Heteroatom. Chem. 9, 651 (1998). 

19. R.N. Barnett, C.L. Cleveland, H. Hakkinen, W.D. 
Luedtke, C. Yannouleas and U. Landman, Eur. Phys. J. 
D9, 95 (1999). 

20. J.M. Soler, M.R. Beltran, K. Michaelian, I.L. Garzon, P. 
Ordejon, D. Sanchez- Portal and E. Artacho, Phys. Rev. 
B61, 5771 (2000). 

21. N.T. Wilson and R.L. Johnston, Eur. Phys. J. D12, 161 
(2000). 

22. T. Li, S. Yin, Y. Ji, G. Wang and J. Zhao, Phys. Lett. 
A267, 403 (2000). 

23. H. Hakkinen and U. Landman, Phys. Rev. B62, 2287 
(2000). 

24. H. Gronbech and W. Andreoni, Chem. Phys. 262, 1 
(2000). 

25. A.H. Romero, A. Castro and A. Rubio, to be published. 

26. M.W. Mathews, R.J. Beuhler, M. Ledbetter and L. Fried- 
man, J. Phys. Chem. 90, 251 (1986). 

27. P.M. Echenique, J.R. Manson and R.H. Ritchie, Phys. Rev. 
Lett. 64, 1413 (1990). 

28. M.H. Shapiro and T.A. Tombrello, Phys. Rev. Lett. 65, 92 



collisions between gold clusters 13 

29. M.H. Shapiro and T.A. Tombrello, Nucl. Inst. Moth. 58, 
161 (1991). 

30. E. Blastein-Barojas and M. R. Zachariah, Phys. Rev. B45, 
4403 (1992). 

31. R. Car and M. Parrinello, Phys. Rev. Lett. 55, 2471 (1985); 
we have used the CPMD code, version 3.5, developed by 
J. Hutter et al, at MPI fur Festkorperforschung and IBM 
Research Laboratory (1990-2001). 

32. J. P. Perdew, K. Burke and M. Ernzerhof, Phys. Rev. Lett. 
175, 531 (1990). 

33. H. Partridge, C.W. Bauschlicher, S.R. Langhoff, Chem. 
Phys. Lett. 175, 531 (1990). 

34. J.M. Seminario, J.M. Tour, Int. J. Quant. Chem. 65, 749 
(1997). 

35. N. Troullier and J. Martins, Phys. Rev. B43, 1993 (1991). 

36. M.S. Daw and M.I. Baskes, Phys. Rev. B29, 6443 (1984). 

37. S.M. Foiles, M.I. Baskes and M.S. Daw, Phys. Rev. B33, 
7983 (1986). 

38. J.N. Murrell, R.E. Mottram, Mol. Phys. 69, 571 (1990). 

39. H. Cox, R.L. Johnston, J.N. Murrell, J. Sol. State Chem. 
145, 517 (1999) and references therein. 



Jose Rogan et al.: Rearrangement collisions between gold clusters 



Cluster 


Method 


Symmetry 


E b [cV] 


R[k] 


v [A 3 ] 




EA 


C2-U 


2.405 


1.937-2.438 


3.317 


3 atoms 


CP 


C2-U 


2.182 


2.607-2.753 


4.225 




WJ 




1.759 


2.780 


4.133 




EA 


D 2d 


2.313 


1.904-2.874 


3.996 




CP (a) 


D 2d 


2.234 


2.570-3.120 


5.857 


4 atoms 


CP (b) 


D 2 h (planar) 


2.517 


2.700 


2.346 




CP (c) 


C2„ (planar) 


2.504 


2.648 


29.759 




WJ 


Td 


2.178 


2.793 


5.004 




EA 




2.685 


2.507 


7.287 




CP (a) 


o h 


2.464 


2.801 


4.672 


5 atoms 












CP (b) 


C 2v (planar) 


2.652 


2.715 


7.706 




WJ 


o h 


2.382 


2.801 


4.290 




EA 


o h 


2.825 


2.536 


5.767 




CP (a) 


o h 


2.523 


2.834 


8.663 


6 atoms 












CP (b) 


D 3h (planar) 


2.870 


2.712 


4.365 




WJ 


0, 


2.574 


2.789 


7.667 




EA 




2.937 


2.577 


5.512 


8 atoms 


CP 


D 2d 


2.857 


2.837 


7.385 




WJ 


D 2d 


2.736 


2.790 


6.792 




EA 


D 4d 


2.999 


2.594 


9.794 


10 atoms 


CP 


D 4d 


2.914 


2.829 


36.225 




WJ 


D 4d 


2.837 


2.778 


30.854 




EA (a) 


lh 


3.027 


2.604 


15.190 




EA (b) 


0§ v 


3.089 


2.693 


21.142 


12 atoms 


CP (a) 


0§ v 


2.879 


2.896 


33.032 




CP (b) 


D 4fc 


2.962 


2.814 


8.780 




WJ 


lh 


2.886 


2.776 


18.395 



